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Large-area  flexible  3D  optical  negative  index 
metamaterial  formed  by  nanotransfer  printing 

Debashis  Chanda1,  Kazuki  Shigeta1,  Sidhartha  Gupta1,  Tyler  Cain1,  Andrew  Carlson1,  Agustin  Mihi1, 
Alfred  J.  Baca3,  Gregory  R.  Bogart4,  Paul  Braun1  and  John  A.  Rogers12* 


Negative-index  metamaterials  (NIMs)  are  engineered  structures 
with  optical  properties  that  cannot  be  obtained  in  naturally 
occurring  materials1-3.  Recent  work  has  demonstrated  that 
focused  ion  beam4  and  layer-by-layer  electron-beam  lithography5 * 
can  be  used  to  pattern  the  necessary  nanoscale  features 
over  small  areas  (hundreds  of  pim2)  for  metamaterials  with 
three-dimensional  layouts  and  interesting  characteristics, 
including  negative-index  behaviour  in  the  optical  regime.  A  key 
challenge  is  in  the  fabrication  of  such  three-dimensional  NIMs 
with  sizes  and  at  throughputs  necessary  for  many  realistic 
applications  (including  lenses,  resonators  and  other  photonic 
components6-8).  We  report  a  simple  printing  approach  capable 
of  forming  large-area,  high-quality  NIMs  with  three-dimensional, 
multilayer  formats.  Here,  a  silicon  wafer  with  deep,  nanoscale 
patterns  of  surface  relief  serves  as  a  reusable  stamp. 
Blanket  deposition  of  alternating  layers  of  silver  and  magnesium 
fluoride  onto  such  a  stamp  represents  a  process  for  'inking'  it 
with  thick,  multilayer  assemblies.  Transfer  printing  this  ink 
material  onto  rigid  or  flexible  substrates  completes  the  fabrica¬ 
tion  in  a  high-throughput  manner.  Experimental  measurements 
and  simulation  results  show  that  macroscale,  three-dimensional 
NIMs  (>75  cm2)  nano-manufactured  in  this  way  exhibit  a 
strong,  negative  index  of  refraction  in  the  near-infrared  spectral 
range,  with  excellent  figures  of  merit. 

NIMs  have  values  of  permittivity  and  permeability  that  are 
simultaneously  negative  over  certain  ranges  of  frequencies1-3. 
Building  on  original  work  with  structures  designed  to  operate  in 
the  microwave  regime,  recent  efforts  have  demonstrated  NIMs  for 
terahertz9,10  and  even  optical11-13  frequencies.  Trilayer  stacks  of 
metal-dielectric-metal  with  total  thicknesses  of  ~100nm,  pat¬ 
terned  using  conventional  planar  techniques  into  two-dimensional 
(2D)  mesh  (fishnet)  layouts  with  characteristic  sizes  in  the  range  of 
100-500  nm  yield  some  of  the  best  examples  of  optical  NIMs  for  the 
NIR  wavelength  range11,14,15.  Recent  theoretical14,16  and  experimen¬ 
tal  work4  suggests  that  stacking  up  many  layers  of  such  2D  NIMs 
produces  cascading  effects  and  a  strong  mutual  inductance  coupling 
between  neighbouring  layers  (LC  resonators).  The  result  is  a 
three-dimensional  (3D)  NIM  with  a  broadband  negative  index  of 
refraction  and  relatively  low  loss17,18.  The  only  reported  structure 
of  this  type  was  fabricated  in  serial  processes  using  focused  ion 
beam  (FIB)  lithography4  to  machine  arrays  of  holes  in  a  uniform 
metal/dielectric  multilayer  coating.  Although  this  method  provides 
the  necessary  submicrometre  resolution  and  excellent  control  over 
out-of-plane  geometries,  its  practical  use  is  restricted  to  overall 
lateral  dimensions  of  less  than  20  pan  (ref.  4).  This  size  limitation, 
taken  together  with  the  complexity  of  the  lithography  tools  and 


their  extremely  slow  patterning  speeds,  make  this  technique 
poorly  suited  to  requirements  for  realistic  applications  including 
superlenses  and  photonic  components,  where  cost,  throughput, 
area  coverage  and  long-range  uniformity  are  important  con¬ 
siderations.  Alternative  approaches  to  3D  nanofabrication,  such  as 
those  based  on  colloidal  self-assembly,  interference  lithography 
and  two-photon  polymerization-based  direct  laser  writing  lack  the 
ability  to  embed  multiple  metal/ dielectric  layers  over  a  3D  space. 
Nanoimprint  lithography19,20  and  various  forms  of  soft  lithogra¬ 
phy21-25,  in  their  standard  embodiments,  offer  the  necessary 
resolution,  but  they  do  not  form  3D  structures  easily.  In  the  follow¬ 
ing,  we  describe  a  3D  transfer  printing  method  that  is  directly  appli¬ 
cable  to  the  fabrication  of  3D  NIMs  with  excellent  optical 
characteristics,  in  ways  that  are  scalable  to  arbitrarily  large  areas 
and  that  are  compatible  with  manufacturing.  We  use  the  resulting 
methods  to  fabricate  3D  NIMs  with  11  layers  and  submicrometre 
unit  cell  dimensions,  over  areas  of  >75  cm2,  corresponding  to 
>105  x  105  unit  cells,  all  with  excellent  uniformity  and  minimal 
defects.  The  area  and  number  of  unit  cells  both  correspond  to 
increases  of  more  than  a  factor  of  2  x  107  over  previous  reports4. 
Multiple  cycles  of  printing  with  a  single  stamp  demonstrate  use 
in  a  manufacturing  mode  at  throughputs  that  are  ~1  x  108  times 
higher  than  those  possible  with  state-of-the-art  FIB  systems 
(~2.5  s  per  unit  cell).  Optical  measurements  and  calculations 
show  negative  index  of  refraction  in  the  near-infrared  (NIR)  spectral 
range,  with  values  as  large  as  R e(n)  ~  -7  at  A  =  2.4  pan  and  high 
figures  of  merit  (FOM)  of  ~8  at  A  =1.95  pan,  indicating 
the  comparatively  low-loss  nature  of  these  large- area  printed 
3D  NIMs.  The  outcomes  compare  favourably  with  small-scale  3D 
NIMs  formed  by  FIB,  where  R e(n)  ^  —1.5  at  A  =  1.775  pan  and 
FOM  «  3  at  A  =  1.8  pan  (ref.  4). 

Figure  1  shows  the  three  main  steps  used  in  the  fabrication:  (i) 
‘inking’  a  stamp  with  a  solid,  multilayer  stack  deposited  by  elec¬ 
tron-beam  evaporation  with  a  collimated,  normally  incident  flux 
of  material,  (ii)  contacting  this  stamp  against  a  target  substrate  to 
‘print’  a  patterned  multilayer  structure  by  transfer  of  material  depos¬ 
ited  on  the  raised  regions,  and  (iii)  removal  of  the  unprinted 
material  from  the  stamp  to  prepare  it  for  another  round  of  printing. 
For  the  cases  described  here,  the  stamps  consist  of  silicon  wafers, 
anisotropically  etched  in  the  required  geometry  through  polymer 
masks  defined  by  soft  nanoimprint  lithography.  Many  stamps  can 
be  produced  with  a  single  imprint  mould,  and  each  stamp  can  be 
used  many  times  for  printing.  The  depth  of  etching  (~1  pm)  in 
the  present  case  exceeds  the  total  thickness  of  the  deposited  multi¬ 
layer  stack  (430  nm),  which  consists  of  alternating  layers  of  silver 
(30  nm)  and  MgF2  (50  nm)  on  top  of  a  sacrificial  layer  of  Si02 
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Figure  1  |  Fabricating  3D  NIMs  by  transfer  printing,  a,  Schematic  of  steps  for  printing,  b,  Top-view  SEM  image  of  a  silicon  stamp  (left;  inset,  magnified 
view),  tilted  view  (52°)  SEM  image  of  a  stack  of  alternating  layers  of  Ag  and  MgF2  on  a  silicon  stamp  (middle;  inset,  magnified  top  view),  cross-sectioned  by 
FIB,  and  a  macroscopic  optical  image  of  a  large  (~2.5  x  2.5  cm)  printed  3D  NIM  (right),  c,  Corresponding  SEM  images  of  a  tilted  silicon  stamp  (left),  an 
eleven-layer  Ag/MgF2  stack  (middle)  and  a  printed  3D  NIM  (right:  inset,  magnified  top  view).  Period  Pof  the  structure  is  850  nm,  and  the  depth-averaged 
widths  of  the  ribs  in  the  fishnet  along  the  x-  and  y-directions  are  635  nm  (Wx)  and  225  nm  (Wy),  respectively.  The  thicknesses  of  the  Ag  and  MgF2  layers 
are  30  and  50  nm,  respectively. 


(50  nm)  designed  to  facilitate  release  during  the  printing.  The 
vertical  dimension  is  defined  by  metal/dielectric  stacks  with  indi¬ 
vidual  thicknesses  that  are  subwavelength  in  dimension.  These 
materials  offer  effective  media  where  the  origin  of  negative-index 
behaviour  emerges  from  capacitive  and  inductive  coupling 
between  the  constituent  layers14,17.  Plan-  and  tilted-view  scanning 
electron  microscope  (SEM)  images  of  a  stamp  before  and  after  mul¬ 
tilayer  deposition  appear  in  the  left  and  middle  frames  of  Fig.  lb,c, 
respectively.  The  inset  in  the  left  frame  of  Fig.  lb  identifies  the  key 
geometrical  parameters  of  the  relief  on  the  stamp.  During  electron- 
beam  evaporation  of  the  multilayer  stack,  deposits  form  not  only 
on  the  flat  raised  and  recessed  regions  of  the  stamps  but  also  on 
protruding  structures  on  the  sidewalls  that  result  from  the  Bosh 
ICP-RIE  etching  process.  These  deposits,  which  appear  most  promi¬ 
nently  near  the  tops  of  the  sidewalls  (Fig.  1,  Supplementary 
Fig.  S3),  are  not  transferred  in  the  printing  process.  The  slightly 
non-vertical  sidewalls  (~6-8°)  in  the  multilayer  structures  on  the 
raised  regions  of  the  stamps  (middle  frames,  Fig.  lb,c)  result  from 
diffusion  and  grain  growth  during  evaporation,  primarily  in  the 
MgF2  layers  (Supplementary  Fig.  S3),  consistent  with  known 
phenomena  in  physical  vapour  deposition26,27  rather  than  from 


incomplete  collimation  of  the  evaporator  systems.  The  sidewall 
angles  are  comparable  to,  or  somewhat  larger  than,  those  achieved 
with  FIB4.  In  the  next  step,  contact  printing  delivers  the  multilayer 
deposits  on  the  raised  regions  of  the  stamp  to  a  target  substrate, 
facilitated  by  removal  of  the  Si02  release  layer  with  dilute  hydro¬ 
fluoric  acid,  as  described  in  the  Methods  and  Supplementary 
Information.  A  printed  structure  (area,  6.5  cm2)  on  a  poly¬ 
dimethyl- siloxane  (PDMS)  substrate  is  shown  in  the  right  frames 
of  Fig.  lb.  This  3D  NIM  includes  11  layers  (Ag  on  top  of  five 
repeats  of  MgF2/Ag)  in  a  fishnet  geometry,  where  the  period  P  is 
850  nm,  and  the  depth-averaged  widths  of  the  ribs  along  the  x- 
and  y-directions  are  635  nm  (Wx)  and  225  nm  ( Wy),  respectively 
(right  frame,  Fig.  lc),  consistent  with  the  geometry  of  the  relief  on 
the  stamp  (left  frame,  Fig.  lb)  and  the  sidewall  angle  of  the  multi¬ 
layer.  Figure  2  provides  a  large-area  SEM  image  (Fig.  2a)  and  mag¬ 
nified  top  (Fig.  2c)  and  cross-sectional  (Fig.  2b)  views  of  regions  of  a 
representative  sample  printed  on  PDMS.  The  printing  process  flips 
the  structure,  as  is  evident  by  comparing  the  cross-sectional  view  in 
Fig.  2b  to  those  in  the  middle  frames  of  Fig.  lb  and  c;  no  other 
changes  are  visible.  This  robust  operation  in  printing  is  owing,  in 
part,  to  the  favourable  mechanical  properties  of  the  multilayer 
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Figure  2  |  Large-area,  printed  3D  NIMs  in  supported  and  free-standing 
configurations,  a,  Large-area  SEM  image  of  a  representative  region  of  a 
printed  3D  NIM.  b,c,  SEM  images  of  cross-sectional  (FIB  milled)  and  top 
views  of  this  structure,  d,  SEM  image  of  a  flexible  3D  NIM  membrane 
formed  by  release  and  subsequent  deposition  on  a  solid  support. 


fishnet  stacks,  as  revealed  in  the  image  shown  in  Fig.  2d  of  a  quasi 
free-standing  3D  NIMs  membrane,  created  by  release  from  a 
stamp  and  solution  deposition  onto  a  solid  support.  The  absence 
of  cracking  or  tearing  suggests  surprisingly  high  levels  of  physical 
toughness  and  mechanical  flexibility.  Although  such  structures  are 
mechanically  robust,  handling  or  otherwise  manually  manipulating 
them  is  difficult  owing  to  their  exceptionally  low  bending  stiffnesses. 
As  a  result,  practical  applications  require  a  supporting  substrate  of 
some  type.  In  addition  to  facilitating  the  printing  process,  this 
mechanics  (Supplementary  Fig.  S4)  might  also  enable  wrapping  of 
curved  surfaces  for  integration  with  refractive  lenses  or  other  sub¬ 
strates  of  interest28-30,  or  folding  to  increase  the  overall  thickness 
of  the  material.  Mechanical  stacking  of  3D  NIMs  like  those  shown 
here  represents  a  conceivable  strategy  for  achieving  even  larger 
thicknesses,  although  overlay  registration  requirements  of  ~50- 
60  nm  create  some  engineering  challenges.  Although  negative- 
index  behaviour  occurs  only  for  normally  incident  light,  mechanical 
flexibility  can  be  useful  to  (i)  facilitate  installation  of  NIMs  into 
systems  that  make  use  of  planar  geometries,  (ii)  enable  mounting 
of  NIMs  on  curved  substrates  in  which  operation  only  occurs  on 
light  that  is  incident  in  locally  normal  configurations  and/or  (iii) 
impart  NIMs  with  a  level  of  mechanical  robustness  that  cannot  be 
achieved  in  rigid,  brittle  formats.  To  demonstrate  the  large-area  seal- 
ability  of  this  process,  we  built  macroscale  3D  NIMs  with  dimensions 
of  8.7  cm  x  8.7  cm,  limited  only  by  the  sizes  of  the  available  stamps. 
Figure  3a  shows  a  large  silicon  stamp  ‘inked’  with  an  1 1 -layer  stack.  A 
representative  3D  NIM  printed  from  such  a  stamp  onto  a  flexible 
PDMS  substrate  is  shown  in  Fig.  3b.  The  fabrication  process  is 


Figure  3  |  Macroscale,  printed  3D  NIMs  and  demonstration  of  use  in  a  repetitive  'manufacturing'  mode,  a,  Macroscopic  optical  image  of  a  10  cm  x  10  cm 

multilayer  deposit  on  a  large-area  silicon  stamp,  b,  3D  NIM  printed  with  such  a  stamp  onto  a  flexible  substrate,  in  a  single  step,  c,  Tilted  view  (~15°) 
macroscopic  optical  images  of  three  different  3D  NIMs  printed  using  a  single  stamp,  d,  Corresponding  representative  small-area  SEM  views  of  these 
three  samples. 
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Figure  4  |  Experimental  measurements  and  simulation  results  for  transmission/reflection  and  refractive  indices  of  3D  NIMs.  a,b,  Experimental  and  FDTD 
results  for  transmission  (T)  and  reflection  (R)  spectra  of  a  three-layer  NIM  monolayer  (a)  and  an  eleven-layer  3D  NIM  (b).  c,  Transmission  spectra  collected 
from  five  different  locations  across  the  entire  area  of  a  8.7  cm  x  8.7  cm,  eleven-layer  3D  NIM.  d,  Corresponding  retrieved  indices  for  three  and  eleven  layers 
showing  a  negative  index  of  refraction  in  the  NIR  band,  e,  FOM  of  three-  and  eleven-layer  3D  NIM  structures,  f,  The  4.48°  NIM  prism  showing 
negative-phase  propagation  at  A  =  1.95  [im.  In  all  cases,  P=  850  nm,  depth-averaged  Wx  =  635  nm  and  Wy  =  225  nm,  Ag  thickness  =  30  nm,  MgF2 
thickness  =  50  nm,  and  background  refractive  index  ns  =  1.2. 


compatible  with  a  range  of  flexible  and  rigid  substrate  materials.  The 
stamps  can  be  reused  to  produce  many  such  NIMs  structures  in  a 
manufacturing  mode.  The  three  frames  of  Fig.  3c  show  three  copies 
produced  using  a  single  stamp,  in  a  repetitive  application  of  the 
process  outlined  in  Fig.  la;  corresponding  representative  small- area 
SEM  views  of  these  three  samples  are  shown  in  Fig.  3d.  Visual  inspec¬ 
tion  of  the  samples  and  high-resolution  images  collected  from  them 
indicate  excellent  spatial  uniformity  and  low  defect  density.  (The  red 
and  rainbow  colours  observed  in  these  images  are  due  to  diffraction.) 

The  optical  properties  can  be  quantified  by  spatially  resolved 
spectral  measurements  of  transmission  and  reflection.  Figure  4 
shows  experimentally  measured  and  finite -difference  time-domain 


(FDTD)  computed  average  transmission/reflection  values  for  unpo¬ 
larized  light  with  3D  NIMs  that  have  geometries  reported  above 
having  three  (Fig.  4a)  and  eleven  (Fig.  4b)  layers.  Partial  infiltration 
of  openings  in  the  mesh  structure  by  the  adhesive  material  (PDMS) 
used  in  the  printing  process  (Fig.  2b)  increases  the  background 
refractive  index  ns.  This  effect  moves  the  optical  response  to  a 
higher  wavelength  range  compared  with  the  case  with  air  as  the 
background  medium14.  We  used  ns=  1.2  in  the  FDTD  simulation. 
If  needed,  this  redshift  can  be  counteracted  by  suitable  choice  of 
fishnet  widths  (Wx  and  Wy)u  or  printing  on  less  sticky  substrates 
or  using  adhesiveless  transfer  techniques31.  The  FDTD  simulations 
used  averages  of  transmission/reflection  separately  computed  with 
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planewave  sources  of  light  having  transverse  electric  (TE)  and  trans¬ 
verse  magnetic  (TM)  polarizations,  to  compare  with  the  unpolarized 
light  used  in  the  spectrometer.  Changes  in  optical  responses  associ¬ 
ated  with  non-vertical  sidewall  growth  can  be  compensated  by 
appropriate  modifications  to  the  lateral  dimensions  of  the  structures 
on  the  stamps.  We  used  stamps  with  dimensions  Wx  ~  565  nm  and 
Wy  ~  170  nm  for  the  eleven-layer  case,  and  Wx  ~  600  nm  and  Wy  ~ 
200  nm  for  the  three-layer  case,  to  achieve  similar  depth- averaged 
dimensions  with  an  8°  sidewall  (Wx  =  635  nm  and  Wy  — 
225  nm).  Both  three-  and  eleven-layer  samples  show  good  agree¬ 
ment  between  experiment  and  simulation.  The  original  theoretical 
work  on  these  types  of  3D  NIMs  suggests  that  a  minimum  of  two 
unit  cells  (five  layers)  is  required  to  achieve  the  necessary  mutual 
inductance  coupling  between  layers14.  Samples  with  two  different 
multilayer  stacks  (three  and  eleven  layers)  enable  comparison  of 
optical  properties  between  a  monolayer  NIM  (three  layers)  and  a 
3D  NIM  (11  layers).  The  results  reveal  the  scaling  of  properties 
with  numbers  of  layers,  and  also  validate  the  correspondence 
between  theory  and  experiment.  The  high  spatial  uniformity  of 
the  properties  of  even  the  largest  3D  NIM  that  we  fabricated  are  cor¬ 
roborated  by  transmission  measurements  collected  at  multiple 
locations  (Fig.  4c).  These  data  and  simulation  results  can  be  used 
to  retrieve  the  impedance  (z)  and  refractive  index  (n)  of  the 
samples,  as  shown  in  Fig.  4d.  The  retrieval  method  used  effective 
medium  approaches,  as  described  elsewhere32-34.  As  confirmation 
of  the  correct  identification  of  the  branches  of  this  inverse 
problem,  we  note  that  analysis  of  the  two  NIM  thicknesses  (that 
is,  three  and  eleven  layers)  yields  the  same  values  for  z  and  n 
(Fig.  4d).  Figure  4d  shows  that  the  real  part  of  the  refractive  index 
R e(n)  is  negative  between  1.75  pan  and  2.4  pan,  with  a  minimum 
of  approximately  —4.2  and  —7.8  for  the  three-  and  eleven-layer 
samples,  respectively.  The  value  is  positive  for  wavelengths  longer 
than  2.4  pm,  consistent  with  causality  and  passive  medium  require¬ 
ments  that  Re(z)  >  0  and  Im (n)  >  0  over  the  entire  spectral  range 
between  1.2  pan  and  3  pan  (Supplementary  Fig.  S6).  The  transmis¬ 
sion/reflection  measurements  and  retrieved  parameters  closely 
match  theoretical  predictions  for  structures  of  this  type14. 
Figure  4d  suggests  that  as  the  number  of  layers  increases  from 
three  to  eleven,  the  real  part  of  the  refractive  index  displays  a 
large  discontinuity  at  2.4  p,m.  As  noted  previously14,  these  two 
branches  can  be  shown  to  result  from  different  spatial  modes14. 
The  losses  can  be  estimated  through  a  FOM  (— Re(n)/Im(n))  calcu¬ 
lation,  where  FOM  is  inversely  related  to  lm(n)  (that  is,  a  measure  of 
material  loss).  In  other  words,  higher  Im (n)  implies  increased  loss 
and  reduced  FOM.  Figure  4e  shows  the  results,  which  indicate 
that  the  three-  and  eleven-layer  structures  have  maximum  FOM 
values  of  ~3.5  (at  2.2  pan)  and  ~8  (at  1.95  pan),  respectively 
(FOM  is  ~2  at  the  largest  value  of  negative  R e(n)  ~  —7.8  for  the 
eleven-layer  NIM  at  a  wavelength  of  2.4  pm).  This  reduction  in 
loss  with  increasing  numbers  of  layers  is  consistent  with  the  predic¬ 
tion  that  cascading  leads  to  strong  magneto- inductive  coupling 
between  neighbouring  functional  layers  through  mutual  induc¬ 
tances  resulting  in  a  broadband  negative  index  of  refraction 
with  low  loss14,17.  We  verified  the  negative  phase  propagation  by 
using  FDTD  to  simulate  the  passage  of  plane  waves  through  a 
4.48°  NIM  prism.  Figure  4f  presents  a  single  snapshot  of  the  time 
evolution  of  the  real  part  of  the  electric  field  component  (Re(Ex)) 
as  a  function  of  space  (x-z);  a  negatively  tilted  field  phase  front 
emerging  from  the  NIM  prism  is  apparent,  as  expected  for  this 
wavelength  range.  Such  NIM  prisms  can  be  fabricated  by  FIB 
milling4  of  a  printed  NIM,  whereas  angled  deposition  or  ion- 
milling  over/ through  shadowing  structures  represent  other,  more 
scalable,  possibilities. 

The  results  presented  here  suggest  that,  when  implemented  with 
automated  tools,  these  simple  printing  techniques  can  be  used  for 
nanomanufacturing  large-area  3D  NIMs  in  a  realistic  manner. 


Easy  access  to  materials  at  such  scales  (more  than  seven  orders  of 
magnitude  larger  than  previous  reports)  should  accelerate  the  adop¬ 
tion  of  these  and  related  NIMs  technologies.  The  same  basic 
approaches  may  also  create  new  opportunities  for  NIMs  designs 
that  cannot  be  achieved  in  any  other  way.  For  example,  folding  or 
multilayer  printing  could  yield  NIMs  with  macroscopic  thicknesses, 
embedded  gain  media,  cavity  defects  or  other  structures.  Combining 
these  printing  procedures  with  other  large -area  lithography  tech¬ 
niques  also  seems  promising.  The  results  may  allow  NIMs  designers 
to  lift  previous  constraints  associated  with  more  conventional 
lithographic  approaches. 

Methods 

Silicon  stamp  fabrication.  The  silicon  was  etched  using  a  Bosch  process 
(etch/passivation:  cycle  time,  5  s/5  s;  RIE  power,  20  W/0  W,  SF6/C4F8  flow  rate, 
35  s.c.c.m./llO  s.c.c.m.  for  constant  ICP  power  of  600  W;  etch  rate,  1  fxm/80  s) 
with  SF6  gas,  to  a  depth  between  1  and  2  pm,  through  a  mask  formed  by  soft 
nanoimprint  lithography.  The  moulds  for  this  process  consisted  of  elements  of 
PDMS  (Dow  Corning,  Sylgard)  formed  by  casting  and  curing  against  a  large- 
area  pattern  of  photoresist  defined  on  a  silicon  wafer  by  deep  UV  projection¬ 
mode  photolithography.  Detailed  descriptions  appear  in  the 
Supplementary  Information. 

Electron-beam  deposition.  Multilayer  stacks  were  deposited  using  a  Temescal 
(FC-1800)  six-pocket  electron-beam  evaporation  system.  For  optimal  collimation,  a 
small  filament  (~6  mm)  was  used.  The  sample  was  mounted  almost  vertically  to  the 
source,  so  that  the  evaporated  flux  front  remained  as  normal  as  possible  to  the 
surface  of  the  stamp.  For  deposition  of  Si02,  Ag  and  MgF2,  average  chamber 
pressures  of  ~1.5  x  10~6  and  deposition  rates  of  ~0.1-0.15  nm  s_1  were  found  to 
yield  the  best  film  quality.The  line  edge  roughness  (LER)  of  the  stamp  and  the  top 
layer  of  the  eleven-layer  stack  were  ~10  nm  and  ~20  nm,  respectively  (the  LER  was 
estimated  from  analysis  of  high-resolution  SEM  images).  As  discussed  in  the  main 
text,  neither  this  LER  nor  the  angled  sidewall  geometry  had  significant  effects  on 
the  properties. 

Printing  the  multilayer  stack.  The  printing  process  involved  careful  introduction  of 
dilute  HF  (5%  in  deionized  water)  onto  the  stamp  to  eliminate  the  sacrificial 
Si02layer.  A  flat  slab  of  PDMS  or  a  PDMS-coated  substrate  was  then  contacted  to  the 
stamp  and  peeled  off  to  complete  the  transfer.  The  3D  NIM  stack  on  PDMS  could,  if 
desired,  be  printed  onto  another  substrate  (Supplementary  Fig.  S2).  A  key  feature  of 
the  procedures  is  that  they  involve  only  low-temperature  steps,  with  pristine, 
unaltered  materials  for  the  NIMs  structures;  they  avoid  chemical,  thermal  or  other 
forms  of  degradation  that  can  be  introduced  by  lithographic  processing,  deep 
etching  and  other  steps  in  alternative  approaches.  These  aspects  allow  compatibility 
with  low-cost  substrates  (such  as  plastics)  and  they  may  also  contribute  to  the 
excellent  optical  properties  of  the  resulting  3D  NIMs,  as  highlighted  above  and 
described  subsequently. 

Optical  measurements.  Transmission/reflection  spectra  were  collected  using  a 
x  10-  x  10,  0.1  numerical  aperture  objective  on  an  optical  microscope  (Hyperion 
2000)  coupled  to  a  Fourier-transform  infrared  spectrometer  (Vertex  70)  and 
outfitted  with  a  spatial  aperture  with  a  diameter  of  3.75  mm.  Reflection  spectra  were 
normalized  to  a  silver  mirror  with  96%  reflectivity.  Transmission  spectra  were 
normalized  using  a  bare  substrate. 

FDTD  simulation.  Transmission  and  reflection  spectra  were  calculated  using 
experimental  parameters  for  the  printed  3D  NIMs  structures,  with  commercial 
FDTD  software  package  (Lumerical,  Lumerical  Solutions  Inc.).  A  Drude  model  was 
used  for  the  dielectric  parameters  of  silver  in  the  FDTD  simulation,  with  a  plasma 
frequency  of  9.0  eV  and  scattering  frequency  of  0.054  eV.  The  scattering  frequency 
was  increased  by  a  factor  of  three  compared  with  that  of  bulk  silver  to  account  for  the 
additional  surface  scattering  loss35.  The  good  agreement  between  experiment  and 
FDTD  results,  which  assume  ideal  geometries  and  zero  LER,  suggests  that  the 
observed  losses  are  intrinsic  to  the  materials  and  stack  geometries  used  here.  Slight 
variabilities  in  the  dimensions  of  the  structures,  the  properties  of  the  constituent 
materials  and  the  levels  of  control  (such  as  uniformity  and  edge  roughness) 
associated  with  the  fabrication  are  the  probable  causes  of  the  10-15%  overall 
discrepancy  between  experimental  observation  and  FDTD  simulation.We  found  that 
transmission  through  the  3D  NIMs  deviates  strongly  (>30%)  from  FDTD 
predictions  for  structures  with  more  than  ~15  layers,  limited  by  the  non-vertical 
sidewall  geometries  discussed  previously.  The  development  of  optimized  physical 
vapour  deposition  techniques  that  can  achieve  more  ideal  sidewall  shapes  is  a  topic 
of  current  work. 
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Supplementary  Information: 

M  ethods 

Silicon  stamp  fabrication.  Diluted  (8%  wt)  photoresist  (SU-8  monomer;  MicroChem)  in 
cyclopentanol  was  spin  cast  (500  rpm  for  5  s  followed  by  2000  rpm  for  30  s)  on  a  silicon  wafer 
(M  ontco  Silicon)  coated  with  silicon  nitride  (-200  nm  using  PECVD).  Soft  imprint  lithography 
with  molds  of  poly(dimethylsiloxane)  formed  patterns  in  the  SU-8.  The  PDM  S  molds  resulted 
from  the  casting  and  curing  techniques  of  soft  lithography  with  masters  consisting  of  silicon 
wafers  with  patterns  of  photoresist  formed  by  deep  ultraviolet  projection  mode 
photolithography  .The  regions  of  the  silicon  nitride  not  protected  by  the  SU-8  were  removed  by 
reactive  ion  etching  (RIE)  (Plasma-Therm)  using  CF4.  Next,  the  exposed  silicon  was  etched 
using  inductively  coupled  plasma  RIE  (ICP-RIE,  Plasma-Therm)  in  a  Bosch  process 
(etch/passivation:  cycle  time  =  5  sec/5  sec,  RIE  power  =  20W/0W,  SF6/C4F s  flow  rate  =  35 
sccm/20  seem  for  constant  ICP  power  =  600  W.  Etch  rate  lp,m/80sec)  with  SF6  gas,  to  a  depth 
between  1  and  2  pin.  The  stamp  fabrication  was  completed  by  removing  the  SU-8  and  silicon 
nitride  by  dipping  into  a  piranha  solution  (FI 2S04:H 2O2  =  3:1)  for  15  min  and  then  into  diluted 
HF  (25%)  for  5  sec  followed  by  sonication  in  acetone  (5  min).  The  detailed  steps  appear  in  Fig. 
SI. 

Printing  the  multilayer  stack.The  printing  process  involved  careful  introduction  of  dilute  H F 
(5%  in  Dl  water)  onto  the  stamp  through  a  stack  (4-5)  of  porous  paper  cloths  (TX  1109 
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T echniC loth  II,  ITWTexwipe),  as  shown  in  Fig.  S2.  After  ~30  min,  the  cloths  were  removed  one 
after  one  and  a  flat  slab  of  PDM  S  (Dow  Corning,  Sylgard)  was  contacted  to  the  NIM  stack  and 
then  peeled  off  to  complete  the  transfer.  The  3D  N IM  stack  on  PDM  S  can,  if  desired,  be  printed 
onto  another  substrate,  as  shown  in  Fig.  S2.  The  stamp  was  first  cleaned  by  rinsing  with  5%  HF 
(diluted  with  Dl  water),  and  then  Dl  water.  The  procedure  was  then  repeated  with  rinsing  in 
acetone  and  then  Dl  water,  followed  by  drying  with  N2. 

Finite  difference  time  domain  modeling.  The  FDTD  simulations  usedaverages  of 
transmission/reflection  separately  computed  with  plane  wave  sourcesof  TE  and  TM 
polarizations,  to  compare  with  the  unpolarized  light  used  in  the  spectrometer.  These  simulations 
used  precise  geometries  extracted  from  images  such  as  the  one  in  Fig.  2b,  but  without  any  LER. 
These  modeling  results  are  consistent  with  those  that  use  vertical  sidewalls  but  depth-averaged 
dimensions.  See  Fig.  S5. 
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Figure  SI  Nanotransfer  printing  based  3D  NIM  fabrication  scheme.  Schematics 
showing  complete  steps  from  stamp  fabrication,  multi-layer  growth  to  transfer 
printing. 
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Figure  S2  T ransfer  of  multi-layer  stack  to  a  target 
substrate.  Steps  showing  release  of  multi-layer  stack 
from  the  stamp  to  the  target  substrate. 
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Figure  S3  Estimation  of  angular  growth.  The  52°  tilted  SEM  view  of  FIB  milled 
single  -400  nm  layer  growth  of  (a)  Ag  and  (b)  M  gF2  on  the  silicon  stamp. 
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Figure  S4  Free  standing  3D  NIM.  The  picked  up 
free  standing  3D  NIM  where  a  mild  adhesive  surface 
(scotch  tape)  was  used  to  pick  up  the  free  standing  3D 
stack. 
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Figure  S5  FDTD  computed  average T/R  with  and 
without  small  angular  (8°)  growth.  FDTD  predicted 
transmission  and  reflection  spectra  of  an  11-layer  3D 
NIM  with  and  without  small  angular  (8°)  growth 
structures  with  P  =  850  nm,  and  depth-averaged  Wx 
=  635  nm  and  Wy  =  225  nm.  Ag  thickness  =  30  nm, 
MgF2  thickness  =  50  nm  and  background  refractive 
index  ns  =  1.2. 
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3 -Layers 


Figure  S6  Retrieved  index  and  impedance  of  3- 
layer  Fishnet.  Retrieve  indices  showing  the  real 
(Re(n))  and  imaginary  (I m(n))  parts  of  the  refractive 
index  and  the  real  part  of  the  impedance  (Re(z)). 
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